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Nomenclature 

A  =  scaling  constant 

E  =  electron  eneigy,  keV 

=  low  energy  limit  for  integration,  keV 
Eo  =  electron  e-folding  energy,  keV 
E,  =  electron  eneigy,  300  keV 
E,  =  electron  energy,  3000  keV 
F I  =  average  daily  electron  fluence  above  energy  Et 

J  =  electron  flux,  cmVs-sr 
J\  =  integral  electron  flux  above  energy  £,,  cmVs-sr 

Ji  =  integral  electron  flux  above  energy  Ei,  cm^/s-sr 

Kf  =  planetary  magnetic  index 

k  =  power-law  exponent 

R  =  ratio  of  the  integral  electron  flux  above  £,  to  that 
above  Ei 

W  =  weight  matrix  between  neurons  in  two  layers  of 
the  neural  network 

“  connection  strength  between  neuron  i  and  neuron  j 
Ar  =  time  interval,  seconds  in  one  day  =  86,400  s 
iKp  =  daily  sum  of  the  planetary  magnetic  index,  K, 
n  =  solid  angle,  sr 

Introduction 

HE  long-term  variability  of  relativistic  electrons  in  the 
Earth's  radiation  belts  determines  the  variability  of  the 
radiation  dose  to  spacecraft  components  and  contributes  to  the 
variability  of  electron  precipiution  into  the  middle  atmosphere. 
The  variability  of  the  dose  influences  spacecraft  design,  and 
Ae  variability  of  the  electron  precipitation  may  have  climatolog¬ 
ical  consequences. 

Energetic  electrons  in  the  Earth’s  radiation  belts  pose  a  signif¬ 
icant  hazard  to  electronic  components  on  spacecraft.'  The  total 
radiation  dose  to  a  lightly  shielded  component  can  be  orders 
of  magnitude  larger  than  typical  susceptibility  levels.  This 
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necessitates  the  use  of  additional  shielding  material  which  adds 
significantly  to  the  total  weight  of  the  vehicle. 

Models  of  the  electron  environment  are  used  to  calculate  the 
total  radiation  dose  for  specific  orbits.  The  results  of  such 
computations,  together  with  susceptibility  data,  are  used  to 
specify  the  minimum  shielding  required  for  spacecraft  electron¬ 
ics.  A  variety  of  models  have  been  developed  for  this  purpose 
by  the  National  Space  Science  Data  Center.  The  electron  model 
in  current  use  is  AE-8^  This  model  is  essentially  static  and  is 
thus  only  valid  for  conditions  similar  to  those  present  when 
the  measurements  were  originally  made.  There  is  some  provi¬ 
sion  for  variability.  AE-8  actually  consists  of  two  models— one 
for  solar  maximum  and  one  for  solar  minimum.  The  solar 
maximum  model  predicts  higher  fluxes  and  thus  is  the  more 
severe. 

Measurements  of  highly  leladvistic  electrons  petfotmed  from 
1979  to  1984  by  the  Los  Alamos  spectrometer  for  energetic 
electrons  (SEE)  spectrometers  on  geosynchronous  spacecraft 
showed  only  sporadic  increases  in  the  electrons  from  1979  to 
1981  around  solar  maximum  but  frequent  increases  during  the 
declining  phase  of  the  solar  cycle  from  1982  to  1984.’  The 
mrxlels  are  thus  contrary  to  those  measurements.  A  serious 
shortcoming  in  the  effort  to  develop  long-term  models  is  the 
lack  of  a  consistent  long-term  set  of  measurements. 

Thome'-’  has  suggested  that  the  precipitation  of  energetic 
magnetospheric  electrons  might  be  a  mechanism  by  which 
solar-magnetospheric  processes  could  significantly  affect  the 
middle  atmosphere  and  perhaps  have  climatological  effects. 
This  idea  was  extended  by  Baker  et  al.‘  using  long-term  mea¬ 
surements  of  energetic  magnetospheric  electrons  and  studied 
further  using  measurements  of  atmospheric  optical  emissions 
in  concert  with  election  measurements  by  Callis  et  al.’-*  Because 
the  magnetospheric  population  of  energetic  electrons  varies  in 
an  irregular  and  poorly  understood  way  with  solar  activity  and 
the  solar  cycle,  it  is  crucial  to  have  a  reliable,  long  term  database 
describing  these  elections  for  correlation  with  atmospheric 
studies. 

We  have  developed  a  neural  network  to  model  the  temporal 
variations  of  the  relativistic  electron  flux  at  geosynchronous 
orbit.’  The  results  reported  here  from  the  neural  network  model 
provide  the  first  estimate  of  the  long-term  variability  of  these 
energetic  magnetospheric  electrons. 

The  model  lakes  as  inputs  10  consecutive  values  of  the 
daily-summed  planetary  magnetic  index  iKp.  The  ouqrut  is  an 
estimate  of  the  daily-averaged  electron  flux  for  energies  >3 
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MeV  for  the  tenth  day.  This  model  provides  an  accurate  simula¬ 
tion  and  forecasting  tool  for  the  geosynchronous  electron  envi¬ 
ronment.  In  particular,  the  model  can  be  used  to  infer  electron 
fluxes  for  time  periods  during  which  direct  measurements  are 
not  available. 

Kp  indices  from  1932  to  date  derived  at  the  Geophysical 
Institute,  Gottingen  University,  are  available  from  the  National 
Oceanic  and  Atmospheric  Administration  (NOAA),  National 
Geophysical  Data  Center.'"  They  provide  over  50  years  of  input 
for  the  neural  network  model.  We  have  mn  the  model  using 
IKp  calculated  from  these  Kp  indices  to  estimate  the  long¬ 
term  temporal  variability  of  the  relativistic  electron  flux  at 
geosynchronous  orbit.  We  have  also  estimated  the  annual  radia¬ 
tion  dose  due  to  electrons  to  a  lightly  shielded  electronic  compo¬ 
nent  at  synchronous  orbit  over  this  same  time  period.  The  results 
are  described  here. 

Neural  Network  Model 

The  neural  network  used  for  this  application  has  been 
described  by  Koons  and  Goraey."  We  will  only  repeat  here  the 
salient  features  that  will  help  the  reader  to  understand  this 
application.  The  neural  network  consists  of  three  layers  of 
neurons.  The  10  neurons  that  comprise  the  first  layer  are  con¬ 
nected  to  the  input,  consisting  of  the  values  of  IKp  for  10 
consecutive  days.  Day  0  is  defined  to  be  the  day  for  which  the 
electron  flux  is  calculated.  The  second  layer  of  neurons,  often 
referred  to  as  the  hidden  layer,  consists  of  6  neurons.  It  is 
common  to  construct  neural  networks  such  that  the  number  of 
hidden  neurons  is  half  the  sum  of  the  number  of  inputs  and 
outputs.  A  single  neuron  is  connected  to  the  output  which 
represents  the  logarithm  of  the  average  flux  of  electrons  for 
day  0. 

Connections  only  exist  between  any  single  neuron  and  the 
neurons  in  the  previous  layer  of  the  network.  Neurons  within 
a  given  layer  do  not  connect  to  each  other  and  do  not  receive 
inputs  from  subsequent  layers.  For  example,  neurons  in  layer 
1  send  outputs  to  layer  2,  and  neurons  in  layer  2  take  inputs 
from  layer  1  and  send  outputs  to  layer  3.  The  connection 
strengths  between  any  two  layers  constitute  the  elements  of  a 
real-valued  matrix  ( W).  The  elemental  values  represent  the 
connection  strength  or  weight  between  neuron  i  (in  one  layer) 
to  neuron  j  (in  the  next  higher  layer).  The  weight  matrices 
are  modified  by  training  using  actual  data,  and  these  matrices 
ultimately  contain  all  of  the  information  relating  the  input  (liKp) 
to  the  output  (the  logarithm  of  the  electron  flux). 

The  network  was  trained  and  tested  using  electron  data  col¬ 
lected  by  a  spectrometer  for  energetic  electrons  instrument.  The 
SEE  sensor  was  designed  and  built  by  the  Los  Alamos  National 
Laboratory.  For  a  description  of  the  instrument  see  Baker  et 
al.’  This  design  has  flown  aboard  a  number  of  geostationary 
satellites.  An  edited  dau  set  covering  the  period  from  April 
19,  1982  to  June  4,  1988  from  one  spacecraft,  1982-019,  was 
provided  for  our  use.  The  data  set  consists  of  daily  average 
count  rates  with  the  galactic  cosmic-ray  background  removed. 

The  network  was  trained  using  count  rates  from  the  high- 
energy  (>3  MeV)  electron  channel.  The  results  were  converted 
to  flux  using  a  geometric  factor  of  0.08  cm^-sr  and  an  efficiency 
of  0.3  for  the  3-MeV  channel.  The  training  data  set  consisted 
of  62  days  of  data  from  July  1 .  1 984,  to  August  3 1 ,  1 984.  The 
training  interval  was  selected  on  the  basis  of  data  continuity 
and  the  occurrence  of  several  large  discrete  flux  enhancements 
within  the  chosen  interval.  The  probability  distribution  function 
for  the  flux  for  the  training  period  is  qualitatively  similar  to 
the  distribution  function  for  the  observations  for  the  entire 
period  from  1982  to  1988  as  shown  in  Fig.  6  in  Koons  and 
Gomey."  However,  the  training  period  has  some  excess  values 
of  high  flux  consistent  with  its  selection  as  a  period  of  significant 
flux  enhancements.  To  obtain  convergence  in  the  neural  net- 
woric,  the  training  criteria  was  set  at  10%  of  the  complete  range 
of  output,  corresponding  to  ~0.5  for  the  logarithm  of  the  flux 
or,  equivalently,  about  a  factor  of  three.  Training  required  2652 


passes  through  the  62  patterns  in  the  training  set.  The  weight 
matrices  and  threshold  neuron  values  for  the  trained  network 
are  given  in  Koons  and  Gorney." 

The  model  was  tested  by  comparing  the  model  outputs  with 
measured  fluxes  over  the  entire  six-year  period  from  April  19, 
1982  to  June  4,  1 988.  For  the  entire  data  set  the  rms  logarithmic 
error  of  the  neural  network  is  0.76  and  the  average  logarithmic 
error  is  0.58. 

Radiation  Dose  Results 

The  annual  radiation  dose  due  to  electrons  to  a  lightly 
shielded  electronic  component  at  geosynchronous  orbit  is  taken 
to  be  the  total  annual  dose  from  electrons  with  energies  greater 
than  300  keV.  We  can  use  the  results  from  the  neural  network 
model  to  estimate  this  number  under  the  assumption  that  the 
shape  of  the  average  energy  spectrum  above  300  keV  is  indepen¬ 
dent  of  time.  To  obtain  the  annual  dose,  we  first  estimate  the 
dose  for  each  day  for  energies  >3  MeV  from  the  neural  network 
model.  We  then  scale  this  dose  to  energies  >300  keV  and 
finally  sum  over  the  days  in  the  year.  We  have  used  elecMn 
spectral  data  from  the  literature  to  determine  the  proper  scaling. 

Figure  1  shows  observed  average  electron  spectra  from  10 
keV  to  10  MeV.  The  energies  between  the  dotted  vertical  lines 
in  the  figure  are  the  energies  of  interest  to  this  paper.  The  figure 
is  a  composite  of  data  from  the  rapid  scan  particle  detector 
(RSPD)  on  the  SCATHA  satellite,"  the  high  energy  particle 
spectrometer  (HEPS)  on  SCATHA,"  the  charged  particle  ana¬ 
lyzer  (CPA)  on  satellite  1979-053,'^  and  the  spectrometer  for 
energetic  electrons  on  1 979-053."  The  RSPD  data  were  obtained 
from  Fig.  10  of  Ref.  II.  They  are  the  differential  number  flux 
as  a  function  of  energy  for  perpendicular  elecuons  averaged 
over  all  Kp  and  L  shells  for  00-01-h  LT  (squares)  and  12-13- 
h  LT  (triangles).  The  data  are  comprised  of  measurements  from 
118  ^ys  during  the  period  from  Feb.  28,  1979  to  June  14, 
1980.  The  HEPS  data  were  obtained  from  Table  38  and  Fig. 
58  in  Ref.  13.  It  consists  of  the  average  energy  spectrum  at  L 
=  6.75  from  74  days  between  Feb.  1 1, 1979  and  Feb.  16,  1980. 
Although  the  day-to-day  variation  is  very  large,  Mullen  and 
Gussenhoven'"  found  the  correlation  coefficients  for  a  power- 
law  fit  to  the  average  energy  spectra  exceeded  0.99  for  the 
energy  channels  at  0.83,  1.2225,  and  2.01 1  MeV.  The  CPA  and 
SEE  data  were  obtained  from  Fig.  7  in  Ref.  3.  The  three  curves 
are  the  energy  spectra  for  June  10,  12,  and  15,  1980.  These 
days  bracket  a  magnetic  storm  with  Kp  reaching  6+  on  June 
11,  1980.  The  low  fluxes  (stars)  were  observed  on  the  day 
before  the  storm.  The  intermediate  fluxes  (diamonds)  were 
observed  on  the  day  following  the  storm.  The  high  fluxes  (hexa¬ 
gons)  were  observed  four  days  after  the  storm.  The  general 
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Fig.  I  Composite  electron  spectra  from  10  keV  to  10  MeV. 
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consistency  among  the  spectra  in  Fig.  1  suggests  that  a  fit  to 
the  data  in  the  0.3-3  MeV  range  will  produce  a  reasonable 
spectral  shape  for  the  average  election  spectrum. 

Baker  et  al.‘  fit  the  high-energy  component  (>1  MeV)  of 
the  CPA  and  SEE  data  to  an  exponential  spectrum  of  the  form 


and  obtained  an  e-folding  energy  £n  of  0.587  MeV.  Mullen  and 
Gussenhoven'^  fit  the  HEPS  data  to  a  power  law  of  the  form 


^s-sr-keV  j 


(2) 


and  obtained  a  value  for  the  exponent  of  k  =  -3.88. 

We  have  fit  the  data  from  240  keV  to  4.0  MeV  with  both 
an  exponential  and  a  power-law  spectral  shape.  The  better  fit, 
shown  in  Fig.  2,  is  achieved  with  the  power  law  with  k  =  “4.5. 
Again  the  energies  between  the  dotted  vertical  lines  in  the 
figure  are  the  energies  of  interest  to  this  paper.  The  correlation 
coefficient  r  is  0.975.  For  the  exponential  fit  Eo  =  0.31 1  MeV 
and  r  =  0.959.  This  value  for  £o  is  considerably  smaller  than 
the  value  of  0.578  MeV  obtained  by  Baker  et  al.‘  for  the  fit  to 
the  higher  energy  portion  of  the  spectrum. 

The  integral  flux  above  energy  Ei  is  given  by  the  double 
integral  over  solid  angle  and  energy 

J=rf  A£*d£dn  (3) 

Jt 


y  = -4itA£r'/(k-l-  1)  (4) 

The  ratio  R  of  the  integral  flux  above  energy  £|  to  the  integral 
flux  above  energy  £;  is  then  given  by 


£  =  y,/yj  =  (£,/£!)“'  (5) 

If  £,  is  300  keV  and  Ei  is  3  MeV,  then  the  ratio  for  k  = 
-4.5  is  3160.  If  the  instantaneous  integral  flux  above  3  MeV 
Jj  is  known,  then  the  average  daily  fluence  £|  >  300  keV  is 
given  by 

£,  =  RJj  AflAr  (6) 

where  is  the  instantaneous  integral  flux  averaged  over  a  day 
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for  electrons  with  energies  greater  than  3  MeV.  For  Afl  =  4ir 
and  A/  =  86400  s,  f,  =  3.43  X  10*  Ji. 

This  technique  was  used  to  calculate  the  average  daily  fluence 
for  each  day  from  January  1 , 1 933  through  December  3 1 , 1 988. 
The  annual  fluences  were  then  calculated  by  summing  the  daily 
fluences  over  all  of  the  days  in  each  year.  The  resulting  annu^ 
fluences  are  shown  as  a  function  of  time  as  the  solid  curve  in 
Fig.  3.  Figure  4  shows  the  ratio  of  the  annual  fluence  to  its 
average  value  for  the  entire  period  from  1933  through  1988. 
The  values  range  from  0.25  to  3.4  times  the  average  value 
which  is  9.16  X  10”  cm"’.  The  probability  distribution  function 
for  the  annual  fluence  is  shown  in  Fig.  5.  The  total  width  of 
the  distribution  is  about  an  order  of  magnitude,  spanning  the 
range  from  3.5  X  10”  to  3.5  X  10'*  cm 

The  early  spacecraft  measurements  suggested  that  the  annual 
electron  fluences  were  related  to  the  sunspot  cycle.'*”  The 
annual  mean  sunspot  numbers  for  the  years  from  1933  to  1988 
are  shown  as  the  dotted  curve  in  Fig.  3.  A  comparison  of  the 
data  sets  plotted  in  Fig,  3  shows  that  the  temporal  behavior  of 
the  electron  fluences  only  weakly  resembles  that  of  the  sunspot 
cycle.  A  scatter  plot  of  the  login  of  the  annual  electron  fluence 
vs  the  mean  annual  sunspot  number  is  shown  in  Fig.  6.  The 
solid  line  shows  the  best  linear  fit  to  the  data.  It  is  given  by 

log  10  fluence  =  1 3.82  +  0.002  X  SSN  (7) 

where  SSN  is  the  annual  mean  sunspot  number.  The  correlation 


Year 

Fig.  3  Annual  fluence  of  electrons  with  energies  >300  keV  at 
geosynchronous  orbit  as  a  function  of  time  calculated  using  the 
neural  network  model,  solid  curve;  mean  annual  sunspot  number 
as  a  function  of  time,  dotted  curve. 
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Fig.  2  Power-law  fit  to  the  spectra  in  Fig.  1  over  the  energy  range 
from  240  keV  to  4  MeV. 


Year 

Fig.  4  Ratio  of  the  annual  fluence  of  electrons  with  energies  >300 
keV  to  their  average  annual  fluence  for  the  years  of  1933-1988. 
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Fig.  5  Probability  diitributian  function  of  the  annual  fluence  of 
electrons  with  energies  >300  keV  at  geosynchronous  orbit. 


Fig.  7  Comparison  of  the  annual  fluence  of  relativistic  electrons 
with  energies  >3  MeV  as  estimated  by  using  the  neural  network 
model  dashed  curve;  and  as  observed,  solid  curve. 


Mean  Sunspot  Number 

Fig.  6  Scatter  plot  of  the  annual  fluence  of  electrons  with  energies 
>300  keV  vs  mean  annual  sunspot  number. 


coefficient  between  the  logio  flux  and  the  mean  sunspot  number 
is  only  0.35. 

Although  there  is  a  clear  trend  toward  higher  annual  fluxes 
with  increasing  sunspot  number  the  variation  is  very  wide,  and 
one  of  the  highest  predicted  fluences  occurs  at  one  of  the  lowest 
sunspot  numbers. 


used  for  the  comparison  in  Fig.  7.  The  neural  network  estimate 
is  systematically  high  by  a  factor  of  2.4.  The  correlation  coeffi¬ 
cient  between  the  estimated  and  the  observed  fluence  is  0.64. 
This  relatively  low  correlation  is  primarily  caused  by  a  large 
discrepancy  (a  factor  of  5.6)  in  1983.  Excluding  that  year  the 
correlation  coefTicient  is  0.97. 

Discussion 

Accuracy  of  Radiation  Belt  Models 

A  curious  feature  of  the  annual  fluences  predicted  by  the 
model  is  the  long  period  of  low  fluences  shown  in  Fig.  3  in 
1961-1972.  Unfortunately,  most  of  the  observations  that  were 
used  to  produce  the  electron  models  by  NSSDC  were  made  by 
spacecraft  during  this  time  period.  The  low  fluence  values 
shown  in  Fig.  3  during  the  1960s  suggest  that  the  electron 
models  might  be  low  by  about  a  factor  of  two. 

Future  Work 

The  main  limitation  of  the  present  work  is  the  assumption  that 
the  shape  of  the  differential  energy  spectrum  of  the  relativistic 
electrons  is  independent  of  time  and  on  average  is  given  by  a 
power  law  with  k  =  -4.5  in  the  range  from  300  keV  to  3 
MeV.  Further  work  is  required  to  determine  the  validity  of  this 
assumption.  The  flux  is  observed  to  be  highly  variable,  and  the 
spectrum  tends  to  harden  as  the  flux  increases.  This  is  apparent 
in  the  data  shown  in  Fig.  1  above  2  MeV.  It  must  sull  be 
determined  if  there  is  a  long-term  secular  variation  to  the  shape 
of  the  average  energy  spectrum  in  the  energy  range  discussed 
in  this  paper. 


Comparison  with  Measurements 
The  results  from  the  neural  network  model  agree  qualitatively 
with  the  measurements  of  relativistic  electrons  performed  from 
1979  to  1984  by  the  Los  Alamos  SEE  spectrometers  on  geosyn¬ 
chronous  spacecraft.  These  measurements  showed  only  spo¬ 
radic  increases  in  the  electrons  from  1979  to  1981  around  solar 
maximum  but  frequent  increases  during  the  declining  phase  of 
the  solar  cycle  from  1982  to  1984.’  The  neural  network  model 
produces  logm  electron  fluences  of  14.01.  13.63,  and  14.09  for 
1979-1981  and  14.47,  14.27,  and  14.29  for  1982-1984.  Thus 
the  earlier  three  years  were  all  near  the  low  end  of  the  probability 
distribution  function  shown  in  Fig.  5  whereas  the  latter  three 
years  were  near  the  high  end. 

Figure  7  shows  a  comparison  of  the  estimate  of  the  annual 
fluence  from  the  neural  network  model  with  the  observed  aver¬ 
age  annual  fluence  determined  from  the  observations  for  the 
years  1983-1987.  The  unsealed  data  at  energie.s  >3  MeV  were 
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TECHNOLOGY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect-engineer"  for  national  security 
programs,  specializing  in  advanced  military  space  systems.  The  Corporation's  Technology 
Operations  supports  the  effective  and  timely  development  and  operation  of  national  security 
systems  through  scientific  research  and  the  application  of  advanced  technology.  Vital  to  the 
success  of  the  Corporation  is  the  technical  staffs  wide-ranging  expertise  and  its  ability  to  stay 
abreast  of  new  technological  developments  and  program  support  issues  associated  with  rapidly 
evolving  space  systems.  Contributing  capabilities  are  provided  by  these  individual  Technology 
Centers; 

Electronics  Technology  Center:  Microelectronics,  solid-state  device  physics, 

VLSI  reliability,  compound  semiconductors,  radiation  hardening,  data  storage 
technologies,  infrared  detector  devices  and  testing;  electro-optics,  quantum  electronics, 
solid-state  lasers,  optical  propagation  and  communications;  cw  and  pulsed  chemical 
laser  development,  optical  resonators,  beam  control,  atmospheric  propagation,  and 
laser  effects  and  countermeasures;  atomic  frequency  standards,  applied  laser 
spectroscopy,  laser  chemistry,  laser  optoelectronics,  phase  conjugation  and  coherent 
imaging,  solar  cell  physics,  battery  electrochemistry,  battery  testing  and  evaluation. 

Mechanics  and  Materials  Technology  Center:  Evaluation  and 
characterization  of  new  materials:  metals,  alloys,  ceramics,  polymers  and  their 
composites,  and  new  forms  of  carbon;  development  and  analysis  of  thin  films  and 
deposition  techniques;  nondestructive  evaluation,  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  development  and  evaluation  of 
hardened  components;  analysis  and  evaluation  of  materials  at  cryogenic  and  elevated 
temperatures;  launch  vehicle  and  reentry  fluid  mechanics,  heat  transfer  and  flight 
dynamics;  chemical  and  electric  propulsion;  spacecraft  structural  mechanics,  spacecraft 
survivability  and  vulnerability  assessment;  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  lubrication  and 
surface  phenomena. 

Space  and  Environment  Technology  Center:  Magnetospheric,  auroral  and 
cosmic  ray  physics,  wave-particle  interactions,  magnetospheric  plasma  waves; 
atmospheric  and  ionospheric  physics,  density  and  composition  of  the  upper 
atmosphere,  remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared 
astronomy,  infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere;  effects 
of  electromagnetic  and  particulate  radiations  on  space  systems;  space  instrumentation; 
propellant  chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
atmospheric  chemical  reactions,  atmospheric  optics,  light  scattering,  state-specific 
chemical  reactions  and  radiative  signatures  of  missile  plumes,  and  sensor  out-of-field- 
of-view  rejection. 


